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(57) Abstract 

A highly efficient Reed-Solomon decoding algorithm and a VLSI implementation thereof are used to correct errors found in data which 
was previously encoded utilizing Reed-Solomon encoding. The major functional units in the Reed-Solomon decoder arc a syndrome 
computation unit (12), a modified Berlekamp-Massey algorithm unit (14), a polynomial evaluation unit (16), a Chien search unit (18), and a 
delay register (10). The Reed-Solomon decoding takes advantage of certain desirable features of Euclid's algorithm and of the Berlekamp- 
Massey algorithm to provide a simplified and computationally efficient algorithm which eliminates the need to perform inverses. 
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(57) Abstract 



A highly efficient Reed-Solomon decoding algorithm and a VSLI implementation thereof are used to correct errors found in data which 
was previously encoded utilizing Reed-Solomon encoding. The Reed-Solomon decoding algorithm takes advantage of certain desirable 
features of Euclid's algorithm and of the Berlekamp-Massey algorithm to provide a simplified and computationally efficient algorithm 
which eliminates the need to perform inverses. 
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REED-SOLOMON DECODER AND VLSI IMPLEMENTATION THEREOF 

RELATED APPLICATIONS 

This patent application claims the benefit of the filing dates of United States Provisional 
10 Patent Application Serial Number 60/055,275, filed on August 13, 1997 and entitled 
SIMPLIFIED PROCEDURE FOR CORRECTING ERRORS OF REED-SOLOMON CODE 
USING AN INVERSIONLESS THE [SIC] BERLEKAMP-MASSEY ALGORITHM; and 
United States Provisional Patent Application Serial Number 60/072,407, filed on January 23, 
15 1 998 and entitled VLSI DESIGN OF A DECODER FOR CORRECTING ERRORS OF REED- 
SOLOMON CODES USING AN INVERSIONLESS BERLEKAMP-MASSEY ALGORITHM, 
the contents of both of which are hereby incorporated by reference. 

on FIELD OF INVENTION 



The present invention relates generally to forward error correcting methodology. The 
present invention relates more particularly to a method for decoding Reed-Solomon encoded data 
so as to effect error correction in a manner which is computationally efficient and which is suited 
for very large scale integrated circuit (VLSI) implementation. The present invention further 
comprises a VLSI implementation of the Reed-Solomon decoding methodology. 

BACKGROUND OF THE INVENTION 

As data storage densities and data transmission rates increase, the ability of hardware 
devices to correctly recognize binary data diminishes. Such binary data, which comprises a string 
of bits, i.e., zeros and ones, must be correctly recognized so as to facilitate the proper 
interpretation of the information represented thereby. That is, each individual zero and one needs 
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to be interpreted reliably, so as to facilitate the accurate reconstruction of the stored and/or 
transmitted data represented thereby. 

Such increased storage densities and transmission rates place a heavy burden upon the 
hardware utilized to recognize the state of individual bits. In contemporary storage devices, high 
storage densities dictate that the individual data bits be stored very close to one another. Further, 
in some data storage systems, the difference between the manner in which a one and a zero is 
stored in the media is not great enough to make the reading thereof as reliable as desired. 



High data transmission rates imply either that the time allocated for each bit is 
substantially reduced, or that a modulation scheme is utilized wherein small differences in 
amplitude, frequency, and/or phase indicate a distinct sequence of bits, e.g., such as in multi-bit 
modulation. In either instance, the ability of hardware devices to accurately and reliably interpret 
2Q such transmitted data is substantially reduced. 

Although modern data storage systems are generally capable of limiting read errors to 
approximately one bit in ten billion, the extremely high volume of data processing typically being 
25 performed in contemporary data processing systems makes even such a seemingly low error rate 
unacceptable. In many systems, such an error rate would result in at least one error per day. 



It must be appreciated that the incorrect reading of even a single bit of data may have 
highly undesirable effects. For example, the incorrect reading of a single bit of a computer 
program may cause the program to lock up or otherwise function incorrectly. Also, the incorrect 
reading of a single bit of numeric data may cause a number to be misread, which may have 
serious consequences when dealing with financial data, for example. 
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As those skilled on the art will appreciate, increasing data storage density and/or 
increasing data transmission speed inherently increases the bit error rate, i.e., the percentage of 
error bits contained in a message. It is desirable to be able to identify and correct data errors 
which occur during data storage and/or transmission operations. By being able to identify and 
correct such data errors, data may be stored more densely and/or transmitted more rapidly. 

Further, the ability to correct data errors facilitates the use of mechanical devices, e.g., 
record/read heads, which are manufactured according to reduced tolerances and which are 
therefore less expensive. That is, hardware may be utilized which is more prone to introducing 
data errors, but which costs less. 



Moreover, the ability to correct errors facilitates storage and transmission of data in low 
signal to noise environments. Thus, since data errors can be corrected, more noise can be 
tolerated within a storage medium (or rather the data signal can be recorded at a lower level as 
compared to the noise level, thus facilitating higher storage densities). Similarly, more noise can 
also be tolerated in data transmission, thereby facilitating increased data transmission speed. 

The ability to store and transmit data in low signal to noise environments provides 
substantial benefit in satellite/spacecraft communications, compact disc (CD) recording, high 
definition television (HDTV), digital versatile disc (DVD) recording, network communications, 
and the like. In satellite/spacecraft communications, a weak signal (low signal to noise ratio) 
inherently results from the use of low power transmitters necessitated by weight constraints. In 
compact disc (CD) and digital disc (DVD) recording, being able to record in low signal to noise 
environments facilitates increased storage densities. In high definition television (HDTV) and 
network communications, being able to operate in low signal to noise environments facilitates 
communication at increased data transmission rates. 
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Error detection and correction methodologies for determining the existence of errors in 
stored and/or transmitted data and for correcting those errors are well known. Among the 
simplest of these methodologies is the use of data redundancy, wherein the data is repeated 
several times when it is stored or transmitted. When the data is read, a particular data error is 
unlikely to affect more than one occurrence of the data, such that a majority vote results in the 
correct data being recognized. 



For example, a data byte may be stored three times. A data error is likely to affect only 
one of these three stored data bytes. Thus, the two correctly stored data bytes constitute a 
majority with respect to the incorrectly stored data byte, thereby identifying the data error and 
assuring correct interpretation of the stored data. However, as those skilled in the art will 
appreciate, such a majority vote method is extremely inefficient, since the data must be stored 
or transmitted several, e.g., 3, times. It is preferable to utilize an error detection and correction 
20 methodology which requires as few added bits to the message data as possible. This tends to 
minimize storage requirements and maximize transmission speed. 

In an effort to overcome the deficiencies associated with the inefficient methodology of 
25 the majority vote system, various means for providing data correction which add only a few extra 
bits to the data have evolved. For example, the well known use of parity bits provides a much 
more efficient scheme for indicating the presence of a data error, although the use of simple 
parity alone does not facilitate the correction of such errors. When only parity error detection is 
used, then incorrectly transmitted data must be re-transmitted. Thus, although corrupt data is 
readily identifiable, corrupt data which was only present in a storage medium, and thus is not 
available for re-transmission, may not be recoverable. 



Since it does not facilitate error correction without retransmission, parity is not a forward 
35 error correction methodology. In forward error correction methodologies, sufficient information 
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is added to the stored or transmitted data to facilitate the reliable correction of data errors without 
requiring such retransmission, as long as a mathematically determined error rate is not exceeded. 

It is desirable to provide a methodology for correcting burst errors. Burst errors involve 
the corruption of a plurality of consecutive bits. Burst errors are common in data storage and 
data transmission. 

Burst errors typically occur during data storage due to faulty media. For example, a 
scratch or other imperfection upon a magnetic or optical storage medium may result in the loss 
or corruption of many consecutive bits stored upon the medium. 

Burst errors typically occur during data transmission due to noise or static within the 
transmission frequency band. Signal fading and cross-talk may also cause burst errors during 
data transmission. 



The well known Reed-Solomon encoding methodology provides an efficient means of 
error detection which also facilitates forward error correction, wherein a comparatively large 
number of data errors in stored and/or transmitted data can be corrected. 

Reed-Solomon encoding is particularly well suited for correcting burst errors, wherein 
a plurality of consecutive bits become corrupted. For example, the implementation of Reed- 
Solomon encoding currently used to store information upon compact discs (CDs) is capable of 
correcting error bursts containing as many as four thousand consecutive bits. 

Reed-Solomon encoding is based upon the arithmetic of finite fields. A basic definition 
of Reed-Solomon encoding states that encoding is performed by mapping from a vector space 
of M over a finite field K into a vector space of higher dimension over the same field. 
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Essentially, this means that with a given character set, redundant characters are utilized to 
represent the original data in amanner which facilitates reconstruction of the original data when 
a number of the characters have become corrupted. 

This may be better understood by visualizing Reed-Solomon code as specifying a 
polynomial which defines a smooth curve containing a large number of points. The polynomial 
and its associated curve represent the message. That is, points upon the curve are analogous to 
data points. A corrupted bit is analogous to a point that is not on the curve, and therefore is easily 
recognized as bad. It can thus be appreciated that a large number of such bad points may be 
present in such a curve without preventing accurate reconstruction of the proper curve (that is, 
the desired message). Of course, for this to be true, the curve must be denned by a larger 
number of points than are mathematically required to define the curve, so as to provide the 
necessary redundancy. 

IfN is the number of elements in the character set of the Reed-Solomon code, then the 
Reed-Solomon encode is capable of correcting a maximum oft errors, as long as the message 
length is less than N-2t. 



Although the use of Reed-Solomon encoding provides substantial benefits by enhancing 
the reliability with which data may be stored or transmitted, the use of Reed-Solomon encoding 
according to contemporary practice possesses inherent deficiencies. These deficiencies are 
associated with the decoding of Reed-Solomon encoded data. It should be noted that in many 
applications, the encoding process is less critical than the decoding process. For example, since 
CDs are encoded only once, during the mastering process, and subsequently decoded many times 
by users, the encoding process can use slower, more complex and expensive technology. 
However, decoding must be performed rapidly to facilitate real-time use of the data, e.g., music, 
programs, etc., and must use devices which are cost competitive in the consumer electronics 
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market. Thus, it is desirable to provide a VLSI device for performing Reed-Solomon decoding, 
since such VLSI devices tend to be comparatively fast and inexpensive. 

Decoding Reed-Solomon encoded data involves the calculation of an error locator 
polynomial and an error evaluator polynomial. One of two methods is typically utilized to 
decode Reed-Solomon encoded data. The Berlekamp-Massey algorithm is the best known 
technique for finding the error locator. The Berlekamp-Massey algorithm is capable of being 
implemented in a VLSI chip. However, the Berlekamp-Massey algorithm utilizes an inversion 
process which undesirably increases the complexity of such a VLSI design and which is 
computationally intensive, thereby undesirably decreasing the speed of decoding process. 

Also well known is the use of Euclid's algorithm for the decoding of Reed-Solomon 
encoded data. However, Euclid's algorithm utilizes a polynomial division methodology which 
20 is computationally intensive and which is difficult to implement in VLSI. 

In view of the foregoing, it is desirable to provide a method for decoding Reed-Solomon 
encoded data which is easily and economically implemented as a VLSI chip and which is 
^ computationally efficient so as to substantially enhance decoding speed. 

SUMMARY OF THE INVENTION 

The present invention specifically addresses and alleviates the above mentioned 

30 

deficiencies associated with the prior art. More particularly, the present invention comprises a 
highly efficient Reed-Solomon decoding algorithm and a VLSI implementation thereof. The 
VLSI implementation of the algorithm is used to correct errors found in data which was encoded 
utilizing Reed-Solomon encoding and is particularly well suited for low cost, high speed 
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applications. The algorithm of the present invention takes advantage of certain desirable features 
of Euclid's algorithm and of the Berlekamp-Massey algorithm. 

• According to the present invention, the error locator polynomial and the error evaluator 
polynomial are calculated simultaneously and there is no need to perform polynomial, divisions 
within the reiteration process. Further, the need to compute the discrepancies and to perform the 
field element inversions, as required by the prior art, is eliminated. Thus, the complexity of the 
VLSI implementation is substantially reduced and the processing time is shortened. A software 
simulation of the decoding algorithm of the present invention shows that the decoding process 
is faster than that of the original Euclid's algorithm and faster than the Berlekamp-Massey 
algorithm. 

As discussed above, two different methods, Euclid's algorithms and the Berlekamp- 
Massey algorithm, for decoding Reed-Solomon encoded data are well known. Euclid's algorithm 
offers the advantage of facilitating the simultaneous calculation of the error evaluator polynomial 
along with the error locator polynomial. However, Euclid's algorithm is less efficient in practice 
than the Berlekamp-Massey algorithm because Euclid's algorithm performs polynomial divisions 
within the iteration process thereof. 

The Berlekamp-Massey algorithm is the most widely used technique for determining the 
error locator polynomial and error evaluator polynomial. Although the Berlekamp-Massey 
algorithm can be implemented in a VLSI chip, such implementation is made undesirably complex 
by the Berlekamp-Massey algorithm's iteration process and field element inversions. Thus, the 
error locator polynomial and the error evaluator polynomial can be calculated simultaneously 
when using the Berlekamp-Massey algorithm with inverses. 
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According to the present invention, desirable features of Euclid's algorithm and the 
Berlekamp-Massey algorithm are utilized to reduce the complexity of a VLSI implementation 
of the decoding algorithm and to reduce the processing time required for decoding. More 
particularly, the present invention utilizes an iteration process similar to that of Euclid's 
algorithm. According to the present invention, field element inversion computations are 
eliminated as well. According to the present invention, the error evaluator polynomial is 
computed simultaneously with the error locator polynomial. Further simplifying computation, 
polynomial swaps are avoided. Thus, the decoding algorithm of the present invention is 
substantially easier to implement in a VLSI design, as compared to contemporary algorithms such 
1 5 as Euclid's algorithm and the Berlekamp-Massey algorithm. 

The method for correcting errors in Reed-Solomon encoded data of the present invention 
can be practiced using either a general purpose computer or dedicated circuits, e.g., a VLSI 
20 implementation. As used herein, the term dedicated circuit refers to a circuit which facilitates 
Reed-Solomon decoding and which is not used in general purpose computing. 



25 



These, as well as other advantages of the present invention will be more apparent from 
the following description and drawings. It is understood that changes in the specific structure 
shown and described may be made within the scope of the claims without departing from the 
spirit of the invention. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is an overall block diagram of a pipeline (255, 223) Reed-Solomon domain 
decoder according to the present invention; 

Figure 2 is a systolic array for computing the modified Berlekamp-Massey algorithm 

according to the present invention; 

Figure 3 is a block diagram of a Chien search unit according to the present invention; and 
Figure 4 is a block diagram of the iteration circuit and the polynomial calculation circuit 

according to a VLSI implementation of the present invention. 



DESCRIPTION OF THE PREFERRED EMBODIMENT 



The detailed description set forth below in connection with the appended drawings is 
intended as a description of the presently preferred embodiment of the invention, and is not 
intended to represent the only form in which the present invention may be constructed or utilized. 
The detailed description sets forth the functions of the invention, as well as the sequence of steps 
for constructing and operating the invention in connection with the illustrated embodiment. It 
is to be understood, however, that the same or equivalent functions may be accomplished by 
different embodiments that are also intended to be encompassed within the spirit and scope of 
the invention. 



The present invention comprises a decoding method for correcting errors which occur 
in Reed-Solomon (RS) encoded data by using a modified, i.e., inversionless, Berlekamp-Massey 
algorithm. The modified Berlekamp-Massey algorithm of the present invention provides a 
method which is readily incorporated into a VLSI design and which is comparatively fast as 
compared to contemporary Reed-Solomon decoding algorithms. According to the present 
invention, the error locator polynomial and the error evaluator polynomial are obtained 
simultaneously. As a consequence, the decoder can be made more modular, simpler, and is 
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suitable for both VLSI and software implementation. An exemplary software simulation result 
and an example illustrating the pipeline and systolic array aspects of this decoder structure is 
given for a (225,223) Reed-Solomon code. 

According to the present invention, the modified Berlekamp-Massey algorithm is 
11 improved such that it finds the error locator polynomial and the error evaluator polynomial 
simultaneously. In other words, the error locator polynomial and the error evaluator polynomial 
are computed directly by using the modified Berlekamp-Massey algorithm with no inverses. 
Thus, according to the present invention, the computation of inverse field elements is completely 
15 avoided in Berlekamp's key equation. 

The modified Berlekamp-Massey algorithm is used to reduce the hardware complexity 
for the time domain decoder of Reed-Solomon codes of the present invention. As a consequence, 
2Q an efficient VLSI implementation of time domain Reed-Solomon decoder is provided. 

Time Domain Decoder for Reed-Solomon Codes 

The algorithm for decoding Reed-Solomon encoded data according to the present 
^ invention is provided below. As discussed above, this algorithm avoids the need to perform 
inverses (as required by contemporary decoding algorithms). According to the present invention, 
the error evaluator polynomial and the error locator polynomial are computed simultaneously. 
This simplification of the decoding process makes the algorithm of the present invention suited 
for implementation in very large scale integrated circuit (VLSI) design. 

30 

The conventional decoder methodology is now explained in more detail. Let n be the 
block length of an (*, k) RS code in GF( 2 m ) , where k is the number of m-bit message symbols. 
Also, let d be the minimum distance of the code, where d-1 denotes the number of parity 
35 symbols. Thus, k = n-(d-l). Finally, let the code vector be c = (c 0 ,c t ,c 2 ...,c n-l ) . Also, let the 
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error vector be e = (e 0 ,e„e 2 ...,e B _,). The received vector at the input of the RS decoder is 

r = c + e_( ro)/ . >rj ... t r mM ) .Supposethat v errors occur in the received vector r and assume that 

2v £ d - 1 . Then, the maximum number of errors in RS codes which can be corrected is 

f = L(y-l)/2j (1) 

where I* J denotes the greatest integer less than or equal to x . 
It is well known that the syndromes are defined by 

/•0 /.0 i=0 



where « is a primitive element in GF(2 m ) and 5* , is known for 1 < k £ d - 1 . 
But c 0 + c x a 1 + c 2 a 21 ... + e„_ t a tm - tu =0, for 1 £ k < d - 1 .Then, (2) becomes 

S.=±YX (3) 



1*1 



where Jf. and y, are the/ 1 * error amplitude and the /* error location, respectively. 
Define the syndrome polynomial as 



A substitution of (3) into (4) yields 



But 



*.i 1 + X,x 

Thus (5) becomes 
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Now, let the error locator polynomial be defined by 

*(*) = 1^(1 + *,*) = (7) 
10 >=' /=o 

where <r(jt)is the polynomial with zeros at the inverse error locations with <r 0 =l and 
deg{a(*)}=v. 

Finally, let the error evaluator polynomial be defined by 
P(x) = i,Y l X l xYl(l + X t x) 



15 



30 



(8) 



where deg{ P(x) }=* and P(0) = 0 . Hence (6) becomes 

20 S<»).*W mod x" W 

<r(*) 

Thus, by (9), Berlekamp's key equation is given by 

(l + S(x))o-(x) = fi(x) mod x rf (10) 
25 where Q(x) = <r(x) + P(x) with Q(0) = 1 and deg{ ft(x) }= v . 
Thus, the error evaluator polynomial is given by 

/>(x) S Q(x) + <r(x). (11) 



35 



It is well known that the Berlekamp-Massey algorithm can be used to find <y(x) and 
Cl(x) in (10). For a given these known polynomials, the error evaluator polynomial can be 
obtained by (1 1). Note that P{X: 1 ) = Cl(X; l ) + <r(X; l )=Cl(X; 1 ) , where x? is the reciprocal 
root of a(x) for 15/ < v. 
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The following algorithm is the iterative procedure developed by Berlekamp and Massey 
to find a(x) and Cl(x) that satisfies (10). 

1 . " Initially define C<°> ( x ) = 1, D <°> (x) = 1, A {0 \x) =1,5 (0 >(x) = 0, t (0) = 1 and k = 0. 



2. Set k = k + 1 if £ 4 is unknown stop. Otherwise define 

,=o 



(12) 



and let 



C (i >(x) = C (i - , >(x)-A (i > J 4 <t - ,) (x).x 
2><*> (x) = D (M > (x) - A ( *> J B ( *- ,) (x) • x 



^(x) = 



xv^Cc), if A<*> = 0 or if 2^-')>Jfc-l 
C<*- l >(;cyA<*\ if A<*>#0 and 2C<*-»> s *-l 



W I D<*-»)(x)/A» if 



A ( *> = 0 or if 2« ( *-'> > Jfc-l 
A<*>#0 and 2« ( *-'> s A-l 



(13) 
(14) 



(15) 
(16) 

(17) 
(18) 



and 



« (t) (x) = 



{(*-0 if A W =0 or if 2*<*- | >>Jfc-l 
k-p- l \x), if A<*>*0 and 2{<*"»> s *-l 



(19) 
(20) 



3. Return to step 2. 
In this iterative procedure, 

<r(x) = C iJ - % \x) 

and 

Q(x) = D (J -'>(x)- 



(21) 
(22) 
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The errors are located by finding the reciprocal root of the polynomial a(x) in (21). 
These roots are usually found by a Chien-search procedure. Note that in (16) and (18) there are 
only two places in the algorithm where the inversion of the element a (A) is required. It is well 
known that inversion takes many operations or special hardware. The modified Berlekamp 
Massey algorithm is developed for finding the error locator polynomial and the error evaluator 
polynomial that satisfies (10). The following algorithm is the modified Berlekamp-Massey 
algorithm which eliminates the requirement for inversion. 

1. Initially define C (0) (*) = 1,D<°> (*)==!, J (0) (x) = 1,5 (0) (jc) = 0,7 (0) =0, 
r <*>=l if JcZO. 

2. Set k = k + 1 if k = d - 1 , stop. Otherwise, define 

a^Srv, (23) 



and let 



C <*>(x) = r <*-» c< 4 -»oo - a™I<*-'>(x) ■ x 
0 w (x) = r <*-»> D ( '-'>(x) - A™5 <*-'>(*) • x 



5 (i) (x) = 
{ (A) (x) = 



xy^^x), if A w =0 or if 2? < *" ,) >Jfc-l 
C^'^x), if A ( *>*0 and 2? ( *- l) <; *-l 



r& k ~ l \x), ifA w =0orif 2P'" > *-l 
^""(x), if A (4) *0 and 2? ( *-° <; Jfc-1 

if A™=0 or if 2? ( *-'>>it-l 
if A w *0 and 2? ( *-'> s Jfc-1 



y(*) = 



Y ( * _l) , if A (i) =0 or if 2f ( *- ,) >ifc-l 
A ( *>, if A w *0 and 2J ( *" ,) <; A-l 



(24) 
(25) 

(26) 
(27) 

(28) 
(29) 

(30) 
(31) 

(32) 
(33) 



Return to step 2. 
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5 



10 



The following theorem illustrates the principal relations of these two algorithms. 

Theorem L For 0 < n < d - 1 , 

C ( " ) W = ftr (,) C ( ' ,> (jc) (34) 

D l " ) (x) = flr<' ) D"(x) (35) 



A iH \x) = r lm) A M (x) (36) 

15 J ( " ) (^) = r ( » ) 5 < " ) (x) (37) 

and 

= gin) ( 3 g) 

where c<*>(*), D ( *>(*), /*<*><*), s<*>(x) and £<*> are defined in (13)-(20), for n = k. 

20 

Proof. Theorem 1 can be proved by induction. By definition, the relations in (34)-(38) hold for 
n = 0 . Assume (34)-(38) hold for n = * - 1 . For n = k , a substitution of (34) for n = k - 1 into 
(23) yields 

25 _, t , '^1° 4zJ *-* (39) 

where yW* 3 ! for / = 0or-L 
30 A substitution of (34) for n = k - 1, (36) for n = k - 1 and (39) into (24) yields 



A-2 A-2 



c (*) (x) = y (*-»(J] r » C (w> (x)) _ (fJ r w A <*>) . ^^(x)^ 

3 5 ^r^flr^^-^x)-^^^) 



(40) 
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By (13), (40) becomes 

C<*>( X) = r <*->f} r wc (i >(x) = Y\r { °C<»(x) < 4I > 
'—i (.-i 

Thus, (34) is proved. 

Using a technique similar to that used for derivation of (39), once can derive (35). 

Note by (39) that tf*>=0 if A <*> =0 and * 0 if A <*> * 0 because y w * 0 for -1 < i S k -2 . 
Hence, if = a<*> =0 or if 2 =2 *<*-» * - 1 , then, by (26) and (36) for » = * - 1 , one has 

2 (i) (*) = * • A< k -»(x) = jc • r ik - l) A< l -»(x) ■ (42) 
But r w =r M i n (32) md A <»(x)=x-A (k - l \x) in (15). Thus, (42) becomes 
A< k \x)=y<»A< k \x) . If A 1 *' = A<» * 0 and 2 7 ( *-» = 2*<*-'> < /fc - 1 , then, by (27), (34) for 
»» = *-l and (39), one has 

A (k >(x) = C ik - , >(x) = f\y< l >C {k - , >(x) ' < 43 > 



/=-! 



But, by (16) and (33), (43) becomes ^ <*> (x) = y ^ A w (x) . Thus (36) is proved. 
Using a technique similar to that used for derivation of (36), one can derived (37). 
If ~& k) =A (k)= 0 or if 27 ( *-" = 2£ { *- ,) > A - 1 , then 7<*> „ 7<*-» m £<*-■> = ,<*>. 
If A (i) = A ( *> * 0 or if 2?<*-'> = 2 <?< i -'> ^ /fc - 1 , then ?*>=*_?*-»> =yfc = . 
Thus, (38) is proved. Hence, the theorem is proved. 
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From (38) and (3 5), the modified error locator polynomial and error evaluator polynomial 
are given by. 



C«~»(x) = ft ■ C«- X \x) = fi a{x) = (Aa„/kr l9 ...fi*„) (44) 

and 

10 

Z) (rf - , >(x) = > ?.£)<''- , >(^) = /?.Q W = (^.n 0 , > 9.n i ,..^.fi < ,. 1 ) (45) 



^ where 



(46) 



/? = Y\r U) = C |,w, (0) = D <J -'>(0) 
is a field element in GF(2 m ). 

20 

Therefore, by (44) and (45), the error locator polynomial and the error evaluator polynomial are 
obtained by 

«*-22<* (47) 

25 " 
and 

p f v — . > v , , g^(x) , . (48) 
/>(x) = Q(x) - <r(x) = -i-i - <t(x) 

P 

u The errors are located by finding the reciprocal root of the polynomialo-(x) in (47) or c (</ " ,) (jc) 
in (44). The Forney algorithm can be used to find the error magnitudes fromP(x) given in (8). 
This algorithm is illustrated as follows: 

35 
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By [14], evaluating (8) and (1 1) at x] x yields 

P( X? ) = ) = Y, n (1 - X e Xj l ) • (49) 
t*j 

By (47) and (49), the error magnitudes are given by 

Y W> Q(*;') D (J - l) (X-/) (50) 

where a \X~ l ) = XjY[(1 -X t X] 1 ) and C^ 0 '^;') = ,0 -^(jr;*) are the derivative with 
respect to * of a(x) given in(7) and c {d ^\x) given in (44) evaluated at * - xj l , respectively. 

Let us recapitulate the decoding of RS codes for errors using the modified Berlekamp- 
Massey algorithm which eliminates the requirement for inversion. This procedures is composed 
of the following four steps: 

Step 1. If v >r > the best policy is not to decode the message. Otherwise, compute the 

syndrome values S x ,5 2 ,...5^ from (2). If S k = 0 for 1 < k < d - 1 , the received 

word is a codeword and no further processing is needed. 

Step 2. Use an inversionless Berlekamp-Massey algorithm to determine the modified 

error locator polynomial f3 • a(x) and the modified error evahiator polynomial fi d(x) 
from the known S k for 1 < * £ d - 1 . 

Step 3. Compute the roots of j3 ■ a{x) using the Chien search. The roots of J3 • a(x) are 

the inverse locations of the v errors. 



-19- 



WO 99/09694 



PCT/US98/16868 



Step 4. Compute the error values from (50). 

To illustrate the time domain decoder for correcting errors a (15,9) Reed-Solomon triple- 
error correcting code over GF(2 4 ) is now presented, The representation of the filed GF(2 4 ) 
generated by the primitive irreducible polynomial p ( x ) = X A + x + 1 is given in Table 1 as 

follows: 



Table 1 Representations of the elements of Qp(2 A ) generated by a 4 +a-f I = 0 



30 



Zero and 
Power s of a 


Polynomials over Qp ^ 4 ) 


vectors overgp ^ ) 


0 


0 


0000 


a 0 


1 


0001 


a' 


a 


0010 


a 2 


a 2 


0100 


a 3 


a 3 


1000 


a* 


a+1 


0011 


a 5 


a 2 + a 


0110 


a 6 


a +a 


1100 


a 7 


a 3 +a + l 


1011 


a a 


a 2 +1 


0101 


a 9 


a 3 + a 


1010 


a 10 


a 2 + a + 1 


0111 


a 11 


a 3 + a 2 + a 


1110 


a 12 


a 3 + a 2 + a + 1 


1111 


a 13 


a 3 + a 2 +1 


1101 


a 14 


a 3 + 1 


1001 



35 
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10 



15 



For example, consider a (15,9) Reed-Solomon code over gp ^ ) with a minimum distance 

d=7 In this code, v errors under the condition 2 v *6 can be corrected. Thus, the generator 
polynomial of such a (1 5,9) Reed-Solomon code is defined by 



Assume the message symbols are all zeros. The encoded code word, which is a multiple of g(x) , 
is C(x)=0. Written as a vector, the code word is c » (0,0,0,0,0,0,0,0,0,0,0,0,0,0,0) Assume the 
error vector is * = (0,0,a n ,0,0,a 5 ,0,a,0,0,0,0,0,0,0) - Thus, the received vector is r = c + e . By 
(2), the syndromes 5 A for r are 

14 

5 *=g^^ n (^ 2 )*+a 5 (a 5 )*+a(a 7 )* for 1***6' 

20 This yields S x =a n ,S 2 =l,S J SSB a 14 ,S 4 =a 13 , 5 5 =1 and S 6 =# 11 Thus, the syndrome polynomial 
is 

S(x) = a l2 x + jt 2 + a ! V + a ! V + x 5 +a n x 6 . (51) 
In this example, by ( 1 ), the maximum error correcting capability is / = \_{d - 1) / 2 J = 6 / 2 = 3 . 

25 

The modified Berlekamp-Massey algorithm is applied next to S( x ) given in (5 1). By this 
means, the modified error locator polynomial and the modified error evaluator polynomial are 
determined by a use of the modified Berlekamp-Massey algorithm. These are accomplished by 

the recursive formula (24) and (25), given in Table 2 . From Table 2, one can observe that the 

30 _ 

computation terminates at k = 6 and the modified error locator polynomial Q (6) (x) is obtained. 

That is, 

C {6) (x) = a s +a 4 x + ar 2 +aV (52) 
35 where deg{c (6) (x)}=3 and p = c (6) (0) = a 5 . 
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By (47), the error locator polynomial is a (x) = C (6) (x) I p = 1 + a u x + a 1 l x 2 + a 14 * 3 • By the 
derivative with respect to x of in (52), one obtains 

C (6), (*) = a 4 +aV. (53) 

Using the Chien search, the root of C (6) (x) constitute the set { a 1 ,a~ $ .a 1 }• Also from Table 
10 2, the modified error evaluate polynomial is obtained. That is, 

D 6 (x) = (a 5 +a l0 x + aV +aV). (54) 



1 5 Thus, by (50), the error values are 

_ D™{a 2 ) a 5 + a x \a 2 ) + a 5 {a 2 ) 2 + g "(g-*) 3 u 
1 a 2 C w {a 2 ) a 2 (a* + a 4 (a~ 2 ) 2 ) ~ " 

20 

Similarly, one obtains 



D (6) (q- S ) _ « 5 + « ! V- 5 ) + cy 5 (q- 5 ) 2 +Cg "(a- s ) 3 5 
25 2 a- 5 C (6) '(a- 5 ) a" 5 (a 4 +a 4 (a- 5 ) 2 ) " 



and 



30 



3 a- 7 C (6) "(a- 7 ) a- 7 (a 4 +a 4 (a- 7 ) 2 ) " a 



35 
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Table 2 below is an example of the modified Berlekamp-Massey algorithm used to find the 
modified polynomials fi xr(x) and fi &(z) . 





k 






>l (x) 


o"'<,> 


r* i ( x) 








0 


o 


c" fM <z) - 1 




5""'f*» - i 


o | X / ■ U 


r 1 


« ■ 0 


10 


1 






27 1 " - * - 1 

^ (x)»C (x) 

»t 




"TO). . . 


r «•» - r M 

- a» 


7«» . * - r" 


15 


2 


a - Z c ' *w 

/«• 

- Ct S t ♦ Ci 5, 

- 11 ♦ a n a" 


?•« 


7 , '<*>-f.7 0 <x> 




-x*l-x 


r III . 

- a° 


r u - r M 
- 1 




3 


I™ - Xf'^-, 
-a" 


?*« 


A (x) = C (x) 
= a ,5 n 








-2 


20 


4 


r»" » 
♦ * u l - a' 


•nx , (<*"+x)x 


■ «°x ♦ X s 




J U)-xf tr) 
-a tt «»a'x> 


- a a . 


r° - r M 


25 


5 


•*•• • 

« a* 


•^•^♦vx) 

♦a I -(fl ,, jr*x , )i 


2r°.«-i 

^ (x) - C (*) 
« a* ♦ a\r 






-a 2 




30 


6 


a * Z c > *w 

/.a 

n Co + Ci S s 
+ Ci o 4 + Ci Jj 

-a , a ,l +a , l 
+ a M a" + a x .a M 

■ 0 


a* 

«Hr , +a*x+<« , +<iV 


2/ W >*-l 
aa*x+ar , x i 




^(D-xi^W 
-a^x+a^aV 


- a 1 


« 3 
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Progra m Implementation and Simulation Results 

The new Reed-Solomon decoding procedure described in the previous section has been 
verified using a C++ program. The program is implemented to correct v errors for general k) 
Reed-Solomon codes, where n is of the form n = 2 OT - 1 , 1 ^ k < n - 1 , and 0 < v < [(/? - k) 1 2j. 

The error patterns are randomly generated over the positions and magnitudes. An example of 
simulation results for a (255,223) Reed-Solomon code correcting v errors, where v < 16 , is 
given in Table 3. For each v , the computation times are in second, which is averaged from the 
computations of correcting 100 different error patterns. The computation times tend to decrease 
when the number of errors occurred decrease. This is because of the quantity ~g*> in Eqs. (24)- 
(33) is zero. 

Table 3. New decoding method for (255,223) RS code. Each entry is the time in second needed 
for one block computation, which is averaged from 100 computations. 



error 


0 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 


time 


0.17 


0.18 


0.19 


0.20 


0.20 


0.22 


0.23 


0.24 


0.26 


0.28 


0.30 


0.32 


0.34 


0.37 


0.40 


0.43 


0.46 



A VLSI Design of th e Time Domain Reed-Solomon Decoder 

Referring now to Figure 1 , the block diagram of a (255,223) Reed-Solomon time domain 
decoder is shown. The major functional units in this decoder are the syndrome computation 
unit 12, the modified Berlekamp-Massey algorithm unit 14, the polynomial evaluation unit 16, 
and the Chien search unit 18. Also, included is a delay register 10, which facilitates 
synchronization of the decoded message output. 
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The syndrome computation unit 12 accepts the received sequence of received messages 
r 0 > r t r 2U and computes their syndromes. 

The sequence of syndromes S n ,...,S 2 ,5, is shifted serially to the modified Berlekamp- 
Massey algorithm unit 14. The output of this unit is the modified error locator polynomial 
C (t6) (x) = p • a(x) and the modified error evaluator polynomial D m (x) = >0 - fi(jc) » where 
^ = c (l6) (0) . The processing of these polynomials is performed by the Chien search unit 18 and 
the polynomila evaluation unit 16, as discussed below. 

Referring now to Figure 2, the VLSI implementation of this modified Berlekamp-Massey 
algorithm of the present invention is shown. The only input is the sequence of 
syndromess^,...^,^. The maximum order of either polynomial C(x) or D(x) is 

/ = 32/2 = 16 . Thus, one needs only t + 1 = 17 sub-registers to store syndromes for calculating 

A register 30 comprises sub-registers F 0 ,F M ...,F I6 . That is, F= (F 0 ,F { ,...,/;,)■ 

Similarly, the registers C =(C 0 ,C 1V .., C l6 ) and DHD 09 D l9 ...,D l6 ) are used to store the 

coefficients of c (*) and d( x ) given in (24) and (25), respectively. Also, C=(C 0 ,C ly .„,C l6 ) 

and £=(£> 0 ,£> |v ..,£> |6 ) are the temporary storage of register C and Z), respectively. The 

register a =( A 0 >A X Z l6 ) and 5 =( 5* o ? gj ^ ^ ) are used to store the coefficients of a ( x ) and 

B ( x ) given in (28) and (29), respectively. The registers y and L are used to store the values 
of y and J respectively. Also, there is a decision circuit which generates a control signal^ to 
control the switches, and a count K is used to store the updated value of k . The outputs are two 
sequences out of the contents of registers c and d following the operations of this structure. 
The operations of this VLSI structure are described below. 

At the initial state, all switches are at position 1. Assume initially in registers that 
F = 0,C = 0,D = 0,c =1>Z> =1> J = l and 5=0. Also, assume the contents in registers y y L 
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and K are 1, 0 and 0 respectively. During k' h iteration, the contents of counter K is increased by 
one. Next, the new syndrome is shifted into sub-register F, and simultaneously the contents in 
sub-registers F ( are shifted right into subsequent sub-registers F M for i =0,1 ,2,..., 15 so that the 

contents of register F will be timed properly for multiplication with the contents of register c 
for generating a (A) . Simultaneously, the contents of registers C and D are loaded into c and d 
registers, respectively. Then, the upper logic is activated to compute a**> given in (23) and 
simultaneously registers J and b are shifted right by one position to perform x • A { *~ 1) ( x ) and 
X .B (*-»(*) respectively. From the known x . J ik ~ X) and x-ff (M, .thfi lower logic is 
activated to compute c {k) ( x ) and D (k) ( x ) signal stored in registers C and D , respectively. At 
the same time, the controls signal 5 is calculated. If a (,) = 0 or 2 > k - 1 ,i.c , 8 =0, 
switches are still remain at position 1 . Thus, and 7<*> in (32) and (38) can be realized and 
the new updated data of ^<*> and 7<*>, i.e. , ^(*)=^<*-o and 7(*>=7<*-»> are stored in registers 

/and L , respectively. If a U) * 0 and 27 ( *"° < k - 1 , i.e. 5=1, one moves switches to position 
2 so that one replaces the contents of registers a and b by the contents of registers c and d , 
respectively. Also, the contents of registers y and L are replaced by a'*> and k - £ ( *" n > 
respectively. Then, the same procedure is received repeatedly. This sequence of operations stops 
when the content of register K equals d - 1 = 32. Finally, the modified error location polynomial 
and the modified error evaluation polynomial stored in registers c and J) » respectively are 
obtained. 

The register transfer language (RTL) of this inverse-free Berlekamp-Massey algorithm 
is described as follows: 

f 0 : F -0, C -0, D -0, c - 1 > D - 1 , J - 1 , B -0, y - 1 , 1 -0, £ -0 and all switches are at position 
1. 
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T t : k-k+\ if k>d -1 = 32 stop. Otherwise, F 0 -S 0 and simultaneously F M -F, for 
'=1,2,..., 15, C-C = C ( *- ,) ( JC )and £>=])(*-»>(*) 
r 2 :Compute m> ,J-x.I (i -'> and 

1 = 0 

T 2 : C-c (k) (x) = r (k ~ l) C-X k) A ,and D-B^{x)^ r ^ l) B-L ik) D 
T A : If 5 = 0 , go to T x . Otherwise, move all switches to position 2 so that J - c =C a " ,} (x) > 
B~D= ( * )> ^ A 1 ** and I - A: - 1 ( *" 0 . Then, move all switch to position 1 and go to 
Tr 

One output of the modified Berlekamp-Massey algorithm unit, the modified error locator 
polynomial c 16 (jc) = p • <t(jc) is sent to a Chien search unit for generating a binary sequence of 

error locations and computing a sequence of p <j'(a~ 2$s ) a~ 15 \...,p -a\a~ { ) a" 1 
simultaneously. In the Chien search method , by [2] ,one has 

7=0 

7-0 

By (54), if one multiplies the f h term of p . a(a~ {i ~ X) ) by a J for 0 < / < 1 6 , then one obtains the 
corresponding term of p . o-(a" ) • Thus, one only needs a constant multiplier and a register for 
computing each term of p . cr(a~ (M) ) ■ In order to facilitate the implementation of the Forney 
algorithm, p . tr(a"') in (54) can be separated into the odd and even terms. That is, 
35 C (l6) (a-') = > ff-<T(a w ) = A e (a") + A 0 (a- / ) forl£y<255 (55) 
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where 

MO= j?>c7,a-"-%z-' 

>=0.2.- 

is the summation of even coefficients and 

»a_ 

is the summation of odd coefficients. 

The summation of odd coefficients of y?-o-(a")» i.e., A 0 (a - ') in (57) is equal to 
/?-o-(a')ar" ;Thatis » 

A 0 (a"') = P ■ <T\a')a' = C m {a , )a > for 1 * ' * 255 

Referring you to Figure 3, the architecture of the Chien search method unit is shown. As 
shown, the odd and even number terms are summed separately. The architecture of this unit is 
similar to that of a contemporary Chien search unit, except there are only 17 error locator sub- 
cells in this (255,223) Reed-Solomon decoder and a zero detector at the end. The zero detector 
is converted the sequence C (,6) (a-')'j into a sequence of Ps and 0's , where 1 and 0 indicate 

theoccunence,i.e., C (,6) (af') m pa{ar') = 0, or nonoccurrence, i.e., C^V^PoteVO, 

respectively, of and error at a specific inverse location. a -< . 

The other output of the modified Berlekamp-Massey algorithm unit, a sequence of the 
modified error evaluator polynomial 0 Cl(x) is sent to the polynomial evaluation unit to 
perform the evaluation of fi-Cl(x). The output of the polynomial evaluation unit 
y?Q(a"')'jand one output of the Chien search unit , a sequence of 



(56) 



(57) 



-28- 



WO 99/09694 



PCT/US98/16868 



pcf(a^a 25 \..fi-d(a x )a x shifted to a finite field divider for computing the sequence 
of n(a- 255 )/a- 2 " a\a' 2iS )^Cl{a x )la x a'Ca-^.Thepioductofthe Q( a ' i )/a' 1 ^(a'Js 
and the other output of the Chien search unit, the binary sequence of error locations forms the 
sequence of the estimated error patterns, i.e., e^e^e^ . Finally, the estimated code vector is 

obtained by subtracting the error vector e from the received vector r . 

Refering now to Figure 4, a very large scale integrated circuit (VLSI) implementation of 
the present invention is shown. The VLSI implementation comprises an integrated circuit 40 
comprising an iteration circuit 42 which performs iteration according to the modified Berlekamp- 
Massey algorithm of the present invention and a polynomial calculator circuit 44 which defines 
the modified error locator polynomial according to equation 44 and also defines a modified error 
evaluator polynomial according to equation 45. 

It is understood that the exemplary algorithm described herein and shown in the drawings 
represents only a presently preferred embodiment of the invention. Indeed, various modifications 
and additions may be made to such embodiment without departing from the spirit and scope of 
the invention. For example, those skilled in the art will appreciate that the algorithm of the 
present invention is applicable to Reed-Solomon codes having various different sizes. Thus, the 
description herein of a (225, 223) Reed-Solomon code is by way of illustration only, and not by 
way of limitation. Further, those skilled in the art will appreciate that various different VLSI 
implementations of the present invention are contemplated. Thus, these and other modifications 
and additions may be obvious to those skilled in the art and may be implemented to adapt the 
present invention for use in a variety of different applications. 

The following is a C-H- program listing of the computer simulation of the algorithm for 
decoding Reed-Solomon encoded data according to the present invention. 
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//- 

#include <stdio.h> 
#include <stdiib.h> 
//include <string.h> 
#include <iostream.h> 
//include <math.h> 
//include <time.h> 

#define N 256 // array size, use 16 to save time when possible 

"«endl 



#define delimit 



//defineDELIMn" w «endl 



20 



25 



// GF class 

//= 



class GF 

15 ( 

public: 

int power; 

int symbl; 

GF(void); 

GF(int pwr, int sym); 
-GF(void); 

GF operator+(const GF & b) const; 
GF operator*(const GF & b) const; 
GF inv(void) const; 
bool operator=(const GF &b) const; 
bool operatorlKconst GF &b) const; 

}; 



// POLY class 

II 



class POLY 
{ 

public: 
GF c[N]; 
POLY(void); 
30 -POLY(void); 

int deg(void) const; 
POLY operator*(const GF & b) const; 
POLY operatoi4-(const POLY & q) const; 
POLY operator*(const POLY & q) const; 
//POLY operator/(const POLY & q) const; 
POLY operator«(const int & j) const; 
POLY operator»(const int & j) const; 
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bool operator=(const POLY & q) const; 

POLY Ishift(void) const; 
POLY rshift(votd) const; 
POLY diff(void) const; 
POLY mod(const int & a) const; 
GF eva!(const GF & x) const; 
void disp(char* title); • 

}; 



// global variables 

" 

void init(void); 

void gen_GF(void); 

void dump(void); 

void gen_errata(void); 

void Fomey_syndromes(void); 

void syndromes(void); 

void Chien(const POLY & p, int * npt, GF *rtpt); 
void sJf_BM_e(void); // simplified ifBMj* 

void POLYdivide(const POLY & p,consfPOLY & q,POLY & quotient,POLY & remainder); 

POLY pdtmoda(const POLY & p,const POLY & q); 

POLY pdtmod(const POLY & p,const POLY & q,const int a); 

void report(void); 

inline double randomOlQ {return (double)random(32767)/32767.;} 



int n,k,d; // n,k command line parameters, d=n-k+l 

int m; // number of bits, defined in inh() 

int Prim[65]; // 2^, m in [0,64], max 64 bits 

GF GFzero,alpha[N]; 

int s0,t0,b0; 

int t,s,v; 

int Power jab[N],SymblJab[N]; 
int Bits[256]; 

int cnt_s; // used as a counter for the erasure s 

int cnt_t; // used as a counter for image frame or t 

int cnt_b; 
int fsize; 

POLY POLYzero,POLYone; 
POLY g; 

POLY erasure,error,errata; 
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POLY inf ,cw,r,eccw; 

POLY Lambda; // erasure locator poly 

POLY ersjp; // erasure locator poly 

//POLY sigma; 

//POLYtau; 

POLY S; 

POLY t; 

POLY errfound; 

int Nor; 

GF RootsfN]; 

int Noe; 

GF Err!oc[N]; 

int a; 

int Reportjlag; 

clock J ckl,ck2,CKI,CK2; 

clockj ctl jUg,ct2_alg,ct_alg; // collect the time used to compute elp and eep 
clockj ctl_Symd,ca_Symd,ct_Symd; // collect the time used for Symdrome 
clock J ct 1 j:hien,ct2_Chien,ct_Chien; // collect the time used for Chien 
clock J ct l_Mag,ct2 J4ag,ct Jtf ag; // collect the time used to compute Magnitude 

clock J ct l_ert,ct2_ert,ct_ert; // collect the time used to compute Magnitude 
double SNR,q2A; // for modulation 



//= 

// methods for GF 

//= 

GF::GF(void) 
{ 

power=n; 
symbH); 
} 

/ / 

GF::GF(int pwr, int sym) 
{ 

power=pwr; 
symbNsym; 

/ / 

GF::~GF(void) 

{ 

} 

//= 

GF GF::operator+(const GF & b) const 
{ 

GF sum; 
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sum .sy mbl=symbi A b.symbl; 
5 sum.power=PowerJab[sum.symbl]; 
return sum; 

} 

/ / 

GF GF::operator*(const GF & b) const 
{ 

GF pdt; 

ifi(symbl=0 || b.symbl=0) return GFzero; 
u pdtpower=(power+b.power) % n; 
pdtsymbHSymbljabfpdtpower]; 
return pdt; 

} 

// 

boo! GF::operator=(const GF & b) const 
{ 

1 5 if (power=b.power && symbl=b.symbl) return true; 
else return false; 
} 

//= ■ 

bool GF::operator!=(const GF & b) const 
{ 

if (power!=b.power || symbl!=b.symbl) return true; 
else return false; 
20 } 

GF GF::inv(void) const 
{ 

///if(symbl=0) {couK^'Inversion of GFzero!!!";} 
ifi(power=0) return alpha[0]; else return alphafn-power]; 

25 



I I 

II methods for POLY 

II 

POLY::POLY(void) 
{ 

30 int i; 

for(i=0;i<n;i++) c[i]=GFzero; 

> 

POLY::~POLY(void) 

{ 

} 



int POLY::deg(void) const 
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{ 

int i,deg; 
deg=n-I; 
foi{M);i<n;i++) 
{ 

if(c[i].syrnbl==0) deg~; else break; 
} 

return deg; 

} 

// 



POLY POLY::operator+(const POLY & q) const 
{ 

int i; 

POLYtmp; 

1 5 for(N0;i<n;i-M.) tmpx[i]=c[i]+q.c[i]; 
return tmp; 

> 

POLY POLY::operator*(const POLY & q) const 

POLY tmp,sum; 
int i; 

20 

for(iH);i<n;i++) tmp.c[n-l-i]=c[n-l-i]; 
for(i=0;i<=q.deg();i++) 
{ 

sum=sum+tmp* qx[(n- 1 )-i] ; 

tmp=tmp.ishiftO; 

} 

9c return sum; 

} 

POLY POL Y::operator*(const GF & b) const 

POLYtmp; 
int i; 

for(i=0;i<-degO;i++) tmp.cfn- 1 -i]=c[n- l-i]*b; 
30 return tmp; 

> 

booi POLY::operator=(const POLY & q) const 
int i; 

fbr(H);i<fc;H-fr-) 

35 < 

if{c[i]=-q.c[i]) continue; 
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else return false; 

5 } 

return true; 
} 

//= - 

POLY POLY::operator«(const int & j) const 
POLY tmp; 

10 inti; 

for(i=0;i<n-j;i+ +) tmp.c[i]=c[t+j]; // the reset components are FGzero????? 
return tmp; 

} 

/ / 

POLY POLY::operator»(const int & j) const 
POLY tmp; 

15 inti ; 

for(i=0;i<n-j;B-+) tmp.c[n-l-i]=c[n-l-(i+j)]; // the reset components are FGzero 
return tmp; 

} 

//= 

POLY POLY::lshiftO const 
{ 

POLY tmp; 
20 int i; 

for(i=0;i<n-l;i++) tmp.c[i]=c[i+l]; 

tmp.c[n-l]=c[0]; 

return tmp; 

} 

POLY POLY::rshift() const 

25 < 

POLY tmp; 

int i; 

tmp.c[0]=c[n-I]; 
for(i=l;i<n;i-H-) tmp.c[i]==e[i-l]; 
return tmp; 
} 

//— 

30 POLY POLY::diffl[void) const 
{ 

POLY tmp; 
int i; 

for(i=n- 1 ;i>=0;i=i-2) tmpx[i]=GFzero; 
for(i=n-2;i>=0;i==i-2) tmp.c[ij=c[i]; 
tmp=tmp.rshift(); 
return tmp; 
^ ) 
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POLY POLY::mod(const int & a) const 
{ 

POLY tmp; 
int i; 
tmp; . 

for(i=0;i<a;i++) tmpx[n- l-i]*=c[n-l-i]; 

//for(i=0;i<n;i++) 

// { 

// j=(n.lM(n-I-i)%a); 
// tnip.c[j]=tmpx[j]+c[il; 
// } 
return tmp; 

} 

// 

GF POLY::eval(const GF & x) const 
{ 

int i; 
GF sum; 

sum=c[n-l-deg()]; 

for(i=deg(> 1 ;i>=0;i~) sum=sum*x+c[n- 1 -i]; 
return sum; 

} 

/ / 

void POLY::disp(char* title) 
{ 

int i; 

if(n>15) return; 
cout.width(9);cout « title; 
cout«":( rt ; 
for(i=0;i<n-l;i++) 
{ 

if(c[i].power=n) couK^-"; 
else 

{ 

cout.width(2); 
com « c[i].power; 
} 

cout « ", M ; 
} 

if(c[i].power=n) cout«"~"; 
else 

{ 

cout.width(2); 
cout « c[i].power; 
} 

cout«T«endl; 
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/ / 

// init : gen_GF, randomize? 

// 

void init(void) 
{ 

if(n== 8-1) m-3; 
if(n= 16-1) m=4; 
if(n=32-I)m=5; 
if(n= 64-1) m=6; 
if(n==128-l)m=7; 
if(n=256-I)m=8; 
gen_GFO; 
T=POLYzero; 

//randomizeO; // seeding the random numbers 
} 



// generate the Galois field 

// 

void gen_GFO 
{ 

int i; 

int tester; // overflow testing 

Prim[4]=0xi3; // ( 1 001 l)=x A 4+x+l 

Prim[6]=0x43; //( 100 001 i^^+x+l 

Prim[8]=0xl Id; // (1 0001 1 101)=x A 8+x A 4+x A 3+x A 2+l 

tester=l«m; // (1 0000), for m=4, for testing if x A 4 exists 

//int sym_tab[16H 1, X 4, 8, 3, 6,12,1 1, 5,10, 7,14,15,13, 9, 0}; 

int symjab[N]; 

sym_tabf0]=l; 

symjablnjN); 

for(i=l ;i<n;i++) // construct the sym table 
{ 

sym jab[i]=sy m_tab[i- 1 ]« 1 ; 

if(symjab[i] & tester) sym__tab[i] A =Prim[m]; 

} 

for(i=0;i<=n;i++)Symbl_tabfi]=sym^tab[i]; 
for(i=0;i<=n;i++) Power Jab[SymbIjab[i]]=i; 
//for(i=0;i<=n;i-H-) cout«Powerjab[i]« w M ;cout«endl; 

for(i=0;i<=n;i++) alpha[i]=GF(i,Symbljab[i]); 
GFzero=alpha[n]; // zero in the field GF 
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for(i=0;i<=n;i++) POLYzero.c[i]=GFzero; // zero poly 

POLYone=POL Yzero; POL Yone.cfn- 1 -0]=alpha[0]; 

g=POLYzero; 

gx[n-l]=alpha[0]; 

POLY tmp; 

tmp.c[n-2]=alpha[0]; 

for(i=I;i<=d-l;i++) 

{ 

tmp.c[n-l]=alpha[i]; 
g=g*tmp; 

//g.disp("genPOLY M ); 
} 

} 



// errata pattern generator 

//M, 

// erasure locator poly, 

// erasure, error and errata=erasure+error are generated 

/ / 

void gen_errata(void) 
{ 

int ij,loc[N],mag[N]; 
POLY tmp_poly; 

randomizeO; 
// random patters 
// random s and t 
if(tO=-I && s0=-l) 
{ 

s=65535;t=65536; 
while((s+2*t)>(cM)) 
{ 

s=d* rand()/RAND_MAX; 
t=(d/2+ 1 )*randO/RAND_MAX; 

} 

} 

else 
{ 

// looping s and t, random are not used 

s=cnt_s; 

t=cnt_t; 

} 

// gen the locations for errata, no repeated 
loc[0]=n*rand()/RAND_MAX; // max returned is n- 1 
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for(i=i;i<s+t;i++) 
{ 

locti]=n*randO/RAND_MAX; 
for(j=Ou<iO-H-) if(loc0]=!oc[i]) -i; 
} 

// erasure locator poly 
ersJp^POLYzero; ersjp.c[n-l]=alpha[0]; 
tmp_poly=POLYzero; 
for(i=0;i<s;i++) 
{ 

tmp _poly.c[n-l]=alpha[0]; 
tmp _poly .c[n- 1 - 1 ]=alpha[loc[i]]; 
ersJp=ersJp*tmp__poIy; 
} 

Lambda=ersjp; // tmp used erasure loc poty in if_BM_ee 

// generate the magnitudes for errata 
// BYTE-wise magnitude 

for(i=0;i<s+t;i++) mag[i]=n*rand()/RAND_MAX; // gen the magnitudes 
// BIT-wise magnitude a^a*! a A 7 

//for(i=0;i<s+t;i++) mag[i]=m*randO/RAND__MAX; //gen the magnitudes 

// finally, the erasure, error and errata 

erasure=POLYzero; 

error=POLYzero; 

foKi=0;i<s;i++)erasure.ctn-l-Joc[i]]=a!pha[mag[i]]: //the erasure 
for(i==s;i<s+t;i++) error.c[n- i -loc[i]]=a!pha[mag[i]]; // the error 

errata=erasure+error; 

// disp the errata pattern 
if(Report flag) 
{ 

cout« ,, e^asure 5 ="«s«*': ,, ; 
for(i=0;i<s;i++) 

cout« l, ( ,, «!oc[i]« w , ,, «erasure.c[n- 1 -loc[i]].power« w ) M ; 
cout«endl; 

cout« M error ="«t« w : rt ; 
for(i=s;t<s+t;i++) 
cout« t, ("«ioc[i]«", ,, «error.c[n-l-loc[i]].power« w ) "; 
cout«endI; 

} 

} 
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" 

// compute the syndromes 

void syndromes(void) 
{ 

intj; . 
S=POLYzero; 

for(p 1 u<=d- 1 Sx[n- i-j]=r.eval(alpha[j]); 



II compute the Forney syndromes from Lambda and S 

void Fomey_syndromes(void) 
{ 

int ij; 
GF sum; 

T-POLYzero; 
forG=lJ<=d-l-sJ-f+) 
{ 

sunv=GFzero; 

for(i=0;i<=s;i-H-) sum=sum+Lambda.c[(n- 1 )-i]*S.c[(n- l>G+s-i)]; 

T.c[n-l-j]=sum; 

} 

} 



II Chien's search 

// no GFzero roots, no repeated roots 

void Chien(const POLY & p, int *npt, GF *rtpt) 

POLY ptmp,pdiv,pquot,prmnd; 
int i; 

*npt=0; 
Nor=0; 

for(i-0;i<n;i++) Roots[i]=GFzero; 
for(i=0;i<n;i++) 
{ 

if(p.evaI(alpha[i])==GFzero) 
{ 
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Roots[Nor++]=alpha[i]; 
///pdiv==POLYzero; 
///pdivxfn- l]=alpha[i]; 
///pdivx[n-2]=alpha[0]; 
///ptmp=p; 
///while(l) 
///{ 

///POLYdivideCptmp.pdiv.pquo^prmnd); 
///if(prmnd===POLYzero && pquoLeval(alpha[il)===GFzero) 

/// { 

/// cout«"repeated roots "«endl; 

/// Roots[Nor++]=alpha[i]; 
/// ptmp=pquot; 

/// } 

/// else break; 

///} 

} 

} 

//if(pxvai(aIpha[0]).symbl==^)Roots[Nor++]=alpha[0]; 



// the devision of 2 POLY, quotient and remainder are returned 

void POLYdivide(const POLY & p f const POLY & q,POLY & quotient,POLY & remainder) 

POLY tmpq; 
GF leadcoef; 
int deg,degp f degq; 

quotient=POLYzero; 
remainder^; 
degp=p,degO; 
degq=q.degO; 

//cout«"— dividing "«endl; 

//cout«"degp« ,, «degp« M degq J = M «degq«endi; 
if(degp<degq) return; 
leadcoe£=qx[n- 1 -degq]; // make q monic 
//cout« n !eadcoef^ w «Ieadcoef.power«endl; 
tmpq=q* leadcoef.invO; 
//quotient.disp("quotient M ); 
//remainder.disp( M remaind B ); 
for(deg s =degp^legq;deg>=0;deg-) 
{ 

// pick the leading coef as the devision since tmpq is monic 
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quotient.c[n- 1 -deg]=remainder.c[n- 1 -(deg+degq)]; 

//((tmpq<<deg)»remainderx[n- 1 -(deg+degq)l).disp( w shifted M ); 

remainder=remainder+ (tmpq«deg)*remainder.c[n- 1 -(deg+degq)]; 

//quotienLdisp( M quotient M ); 

//remainder.disp( M remaind M ); 

} . 

quotient=quotient*Ieadcoef.invO; // Un-make the monic 
//(quotient*q+remainderfp).disp( ,, divcheck n ); 

//cout«"— divided « M «endl; 

return; 
} 



// the product of 2 POLY mod a 

15 



POLY pdtmoda(const POLY & p,const POLY & q) 

POLY sum; 
int ij; 



sum; 

for(i=0;i<a;i++) 
20 for(j=0y<a-ij++) 

sum.c[n- 1 -(i+j)l=sum.c[n- 1 -(i+j)]+px[n- 1 -i]*q.c[n- 1 -j]; 

return sum; 
} 



//= 



// the product of 2 POLY mod a 

I I 



POLY pdtmod(const POLY & pxonst POLY & q,const int a) 

POLY sum; 
int ij; 

sum=POLYzero; 
for(i : =0;i<a;i++) 

for(j=0u < a-iu ++ ) 

sumx[n- 1 -(i+j)]=sum x[n- 1 -C+j)]+p.c[n- 1 -i]*q.c[n- 1 -j] ; 

return sum; 
} 
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// Compute the error polynominal, output is errfound 

// s_if_BM_e : simplified if_BM_e, locator and evahiator are computed 

lj simultaneously. 

void s_if_BM_e(void) 
{ 

int i j,kk; 

POLY Cb t Db,Ab,Bb; 

POLY CbJ f Db_l,AbJ,Bb 1; // previous data 

POLY Cbj>rirae; 

int Lb,Lb_l; 

GF gammab,gammabj ; 

GF deitab; 

GF beta; 

POLY sigma.Omega; 
GFtmp.tmp inv; 
POLY Y; 

//cout« M WeIcome to s_ifBM_e decoder! ! !"«endl; 

// the syndromes 

syndromesQ; 

//starting the BM 

ctl_a!g=clock(); 

//kk=0; 

Cb=POLYzero; Cb.c[n-l]=aipha[0]; 
Db=POLYzero; Db.c[n-l]=alpha[0]; 
Ab=POLYzero; Ab.c[n-l]=alpha[0]; 
Bb^POLYzero; 
Lb=0; 

gammab=alpha[0]; 
/* 

cout« M k= w «kk« H "« H deltabar= ,# «deltab.power«endI; 

Cb.dispCCbar"); 

Ab.dispCAbar"); 

Db.disp("Dbar w ); 

Bb.disp("Bbar M ); 

cout<< ,, gammaba^= ,, <<ga^unab.power<< ,, "<< w Lbar="<<Lb<<endi<<endl; 
*/ 

for(kk= 1 ;kk<=d- 1 -s;kk++) 
{ 

AbJ=Ab; 
Bb_I=Bb; 
CbJ=Cb; 
DbJ=Db; 
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LbJ=Lb; 

gammab_ I =gammab; 
deitab=GFzero; 

forG=0J<-LbJj4+)de!tab^^^ 

Cb=CbJ ♦gammab+Ab_l.lshiftO*deltab; 
Db=Db_l *gammab+Bb__I.lshiftO*deItab; 

if(deltab=GFzero || 2*Lt»kk-I) 
{ 

Ab=Abj.IshiftO; 
Bb=Bb_MshiftO; 
Lb=LbJ; 

gammab=gammab_l ; 
} 

else 

{ 

Ab=CbJ; 
Bb=Db_l; 
Lb=kk-Lb_l; 
gammab=deltab; 
} 

/* 

cout«"k='»«kk«" "« w deItabar="«deltab.power«endI; 

Cb.dispCCbar"); 

Ab.disp("Abar M ); 

Db.disp("Dbar n ); 

Bb.disp("Bbar"); 

cout« ,, gammabar= ,l «gammab.power« H "« ,, Lbar= ,, «Lb«endl; 

cout«endl; 

*/ 

} 

ct2_alg=clock(); 

(*jilg+=(ct2_alg-ctl_alg); 

// using Chien's search to find the roots of sigma 

Chien(Cb,&Nor,Root$); 

Noe=Nor, 

for(i=0;i<Noe;i++) Errloc[i]=Roots[i].inv(); // inverse locations 
//cout« H Chien , s search : Nor ^^Nor^" : "; 
//for(i=0;i<Nor;i++) cout<<Errloc[i].power<<" ";cout«endi; 

// the error values and the error polynomial Y 
Cbj)rime=Cb.difTO; 
for(i=0;i<Nor;i++) 
{ 

tmp=Errloc[i]; 
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tmp_inv=Err!oc[i],invO; 
Y.c[n-l-tmp.power]= 
Db.eva!(tmpJnv)*(Cb .^rime.evaKtmp^invJ^tmp^invXinvO; 

errfound=Y; 
} 



// report 

void report(void) 
{ 

if(Report_flag) 
{ 

info.disp("Info M ); 

cw.disp("Codeword"); 

erasurcdispC'erasure"); 

erTondispe'error"); 

r.disp( M Received H ); 

Lambda.disp("Lambda"); 

S^spC'Syndrome"); 

T.disp( r, Forney M ); 

cout«endl; 

errata.dispC'Errata"); 

errfound,disp("ErrFound"); 

eccw.disp( ,, ErrC^ted ,, ); 

cout«endI; 

} 

if((cw+eccw)=POLYzero) 
{ 

if(Report_flag) cout«"Error corrected! f! M «endl; 
else 

( 

cout« M Error NOT corrected ***************»********«***********»« en( ij. 

errata.disp( ,, Errata"); 
S.disp("syndrome M ); 
(cw+eccw^dispCwrong?"); 
} 

if(ReportJlag) cout«DELIMIT; 
} 
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// Dump some global data 

void dump(void) 
{ 

int i; 
GF tmp; 

//global parameters 

cout« M n= M «n«" k="«k« H d= M «d«" m= M «m«endl; 

cout<< ,, Primitive~ , <<Prim[m]<<endl; 
// alpha 

for(i=0;i<n+ i cout « aIpha[i].power « M " « alpha[i].symbl « "\n H ; 
cout«endl; 

//addition and multiplication table 
//forG=0J<n+ly++) 



// for(H);i<n+l;i++) 

" { 

U tmp=alpha[i]+alpha[j]; 

/; cout « alpha[i].power « n + " « alpha[j].power « H = " « tmp.power «"\t tt ; 

/; cout « alpha[i].symbl « " + " « alphaOJ.symbl « " = " « tmp.symbl «"\t"; 

// if ((i % 2) =0) cout « endl; 

" } 
//cout «endl « endl; 

//for(j=0J<n+lJ++) 

// for(i=0;i<n+l;t++) 

" I 

M tmp=alpha[i]*alpha[j]; 

11 cout « alpha[i].power « rt * " « alpha[j].power « " - " « tmp.power «"\t"; 

11 cout « alpha[ij.symbl « ■ * " « alpha[j].symbl « " = " « tmp.symbl «"\t M ; 

// if ((i % 2) =0) cout « endl; 

// } 



//cout«endl; 
} 



int main(int argc, char * argv[]) 
{ 

char method[80]; 

int start j>,end_s,start_t,endj,start_b,end_b; 
POLY quotient t remainder; 
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if(argc=l) 

5 { 

cout«"eccn k method s t blocks Rpt M «cndl; 

cout«"ecc 255 239 s_if_BM_e 0 3 100 0"«endl; 

cou^^eport^flag-O, U M «cndI«endI; 
exit(Q); 

} 

n=atoi(argv[l]); 
k=atoi(argv[2]); 
1U d=n-k+i; 

strcpy(method,argv[3]); 
s0=atoi(argv[4]); 

start_s=s0; endjs=s0; 

if(s0=999) {start_s=0; end_s=d- 1 ;} //all 
if(sO— 1) {start_s=0; end_s=0;} //randomized 
t0=atoi(argv[5]); 
15 start_t=tO; end_t=tO; 

if(t0=999) start_t=0; // end_t=(d- l-s)/2; moved to the loops 
if(t0=-l) {start_t=0; endj=0;} 
b0=atoi(argv[6]); 

start_b=l; endj>=b0; 
Report_flag=atoi(argv[7]); 



20 



25 



cout«DELIMIT; 

cout<<"RSC<<n<<V'<<k<<") "«method«" d="«d«" m="«m«endl- 

cout«DELIMIT; 

cout«endI«endI; 

cout«" s t Blocks GenErt Syndnn Algthm Chiens Magntd Tota!T"«endl* 
cout«DELIMIT; 

//= 



CKl=clockO; // global timer 
for(cntjs=start_s;cnt s<=end s;cnt s++) 
{ ~ ~ 

if(t0=999) endjKd-l-cnt_s)/2; 
for(cntj=startj;cnt t<=end tent t++) 
{ ~ " 

ckl=clockO; 

30 foKcntJ>=startJ>;cntJ><=end b;cnt b++) 
{ ~ " 

//== 



// get an info from the array image, which is from a file 

info=POLYzero; 

cw^POLYzero; 

3 ^ POL Ydivide(info,g,quotient^emainder); 
cw=info+remainder; 
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// add the errata pattern 

ctl_ert=clock(); 

gen_errata(); 

ct2__ert=clock0; 

ct_ert+=(ct2_ert-ctl_ert); 

r=cw+errata; 

// RS decoder 

// output is errfound 

s_if_BM_eO; 

// error corrected code word eccw = received + error found 
eccw=r+errfound; 

report(); 

} // count the blocks 

ck2=c!ock0; // block timing 

cout.width(3);cout«s; 

cout.width(3);cout«t; 

cout.width(8);cout«b0; 

coutwidth(8);cout«ct__ert; 

cout.width(8);cout«ct_Symd; 

cout.width(8);cout«cralg; 

cout.width(8);cout<<ct_Chien; 

cout.width(8);cout«ct_Mag; 

cout.width(8);cout«(ck2-ckl)«endI; 

ct_ert=0; 

ct_alg=0; 

ct_Symd=0; 

ct_Chien=0; 

ct_Mag=0; 

if(Report_flag) cout«DELIMIT; 
} // count cnt_t 
} // count cnt_s 

CK2-cIockO; // global timing 
iff Report_flag=0) cout«DELIMTT; 

cout«'Totai time= M ;cout.width(5);cout«(CK2-CK 1 )«end!«endl; 

return 0; 
} 
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CLAIMS: 



1 . A method of using a general puipose computer to correct errors in data stored using Reed- 
Solomon encoding by solving an error locator polynomial and an error evaluator polynomial 
using the step of iterating without performing polynomial division and inverses during iteration. 

2. The method as recited in Claim 1, further comprising the steps of: 

a. using a general purpose computer to define a modified error locator polynomial 
according to the equation: 

C<"-'>(*) = P -C<"-V) = fi. a (x) = (fia 09 pa l9 ^a^ x ); and (44) 

b. using a general purpose computer to define a modified error evaluator polynomial 
according to the equation: 

^V)^^^ (45) 

3. The method as recited in Claim 1, wherein the step of iterating comprises using a general 
purpose computer to calculate the error locator polynomial and the error evaluator polynomial 
simultaneously without computing an inverse in a finite field. 

4. The method as recited in Claim 1 , wherein the step of iterating comprises using a general 
purpose computer to iterate according to a modified Berlekamp-Massey algorithm. 

5. The method as recited in Claim 1, wherein the step of iterating comprises: 

a. using a general purpose computer to initially define: 
C (0) (x) = l,Z) (0) (x) = l,^(x) = 1 ,I<°>(x) = 0, £< 0) =0, 
r <*> = 1 if k £ 0 ; 

b. using a general purpose computer to set k = k + 1 if k = d - 1 , otherwise, defining 

a*) = £?<*-v y ; (23) 

J"0 
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and using a general purpose computer to let 

C <*>(,) = y <*-«> c ( *-'>(x) - A ( * , I (1 -"(x) • x 

7J (i >(*) = r < M > D ( '-'>(*) - A™ J<*-»(jc) • * 14 



(24) 
(25) 



^(*) W= *v4 ( *" l) (x), if A w = 0 or if 2? { *- |) >A-l 
|c (t -'>(x), if A w *0 and 2?<*- ,) s *-l 



(26) 
(27) 



bW( x \= \x-B {k - l \x) t if A w = 0 or if 2? ( *' l) > k-l 
\m- l \x), if A w *0 and a**" s *-l 



(28) 
(29) 



j<*>to- , if A™ = 0 or if 2? ( *- ,) >Jt-l 

W " U-i*" 0 . if A<*>,0 and 2?<*-'> s *-l 



(30) 
(31) 



,(*) = , Y 




if A (A) =0 or if 2{ ( *" ,) > k-l 
if A (i) *0 and 2{ ( *' ,) <; k-l 



(32) 
(33) 



; and 

d. repeating step (a). 
6. The method as recited in Claim 1 , further comprising the steps of : 

a. using a general purpose computer to determine an error locator polynomial 
according to the equation: 



; and 

b. using a general purpose computer to determine an error evaluator polynomial 
according to the equation: 



o-(x) 



P 



(47) 



-a(x) 



(48) 
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7. The method as recited in Claim 1, further comprising the step of using a general purpose 
computer to determine error magnitudes according to the equation: 

Y _ _fl(*7') = _gg)__ D^\X~ l ) (50) 



. J 



wherein <r'(X] l ) = X ,Y{{\- X e X/) and C (J ^'(X: 1 ) = B<r'(X: 1 ) are the derivatives 
with respect to * of a{x) and (pf d '\x) evaluated at * = , respectively. 

8. A method for correcting errors in data stored using Reed-Solomon encoding, the method 
comprising the steps of: 

a. using a general purpose computerto compute syndrome values S i ,S 2 ,„.£,,_, when 

teu using 

c -V.V V « fori <*<</-! 

f"0 fcO f=Q 



unless S k =0 iox\<k<d-\; 

b. using a general purpose computer to determine a modified error locator 
polynomial 0 a(x) and a modified error evaluator polynomial 0 ti(x) from S t for 
l£k£d-l via an inversionless Berlekamp-Massey algorithm; 

c. using a general purpose computer to compute the roots of p ■ cr(x) using a Chien 
search, the roots of fi • cr(x) being inverse locations of the v errors; and 
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using a general purpose computer to compute error values according to the 

n(x; { ) Q(X] 1 ) d«- 1 \x] 1 ) . (50) 

9. A method for correcting errors in data stored using Reed-Solomon encoding, the method 
comprising the step of using a dedicated circuit to iterate without performing polynomial 
division and inverses during iteration, 

1 0. The method as recited in Claim 9, further comprising the steps of: 

a. using a dedicated circuit to define a modified error locator polynomial according 
to the equation: 

C <"-> (x) = p - C<"-'> (*) = fi . a(x) = (J3a 0 , fi* x ) ; and (44) 

b. using a dedicated circuit to define a modified error evaluator polynomial 
according to the equation: 

S ( '-V) = ^Z> ( '-V)^^ (45) 

11. The method as recited in Claim 9, wherein the step of iterating comprises using a 
dedicated circuit to calculate an error locator polynomial and an error evaluator polynomial 
simultaneously without computing an inverse in a finite field. 

12. The method as recited in Claim 9, wherein the step of iterating comprises using a 
dedicated circuit to iterate according to a Beriekamp-Massey algorithm. 

13. The method as recited in Claim 9, wherein the step of iterating comprises: 
a. using a dedicated circuit to initially define: 

C (0) (x) = 1 , D< 0 >(*) = 1 , A m (x) = 1 ,B i0) (x) = 0 ,? 0) = 0 , 
r <*> =1 if k<LQ ; 



d. 

equation: 
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14. 



b. using a dedicated circuit to set k = k + 1 if k = d - 1 , otherwise, defining 



y-o 



and using a dedicated circuit to let 

C < 4 >(x) = C< 4 ->(x)- 3*>7<*-'>(x) -x 

D (i) (x) = / ( *-'> D ( * _l) (x) - ( X ) • x 14 



5 w (x) = 
?<*>(x) = | 



xv^'fcX if A» = 0 or if 2? ( *- ,) >*-l" 

C«- l \x), if A<*>*0 and 2? ( *-'> s *-l 

x-S^'^rX if A»=0 or if 2« <*-'>> A- 1 

D«- l \x\ if A (A) *0 and 2^ l) s k-\ 



if A^O or if 2? ( *- ,) > k-\ 
k-V*- x \ if A»*0 and 2? ( *" ,) s A-l 



y(*) = 



y»-i) t if = or if 2« { *- ,) >A-l 
A (i) , if A w *0 and 2{ ( *-'> 



<. k-l 



(23) 



(24) 
(25) 

(26) 
(27) 

(28) 
(29) 

(30) 
(31) 

(32) 
(33) 



a. 

the equation 



; and 

d. repeating step (a). 

The method as recited in Claim 9, further comprising the steps of : 

using a dedicated circuit to determine an error locator polynomial according to 



o-(x) = - 



"(</-D 



(*) 



(47) 



; and 
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b. using a dedicated circuit to determine an error evaluator polynomial according to 
the equation: 

1Q 15. The method as recited in Claim 9, further comprising the step of using a dedicated circuit 
to determine error magnitudes according to the equation: 

Y . = —J ' - '^~i> - (50) 



J no-^;') ^; , o• , (^; l )~^; l c (, '- , >•(x; , ) 



15 

wherein a'(Xj l ) = J^n* 1 311(1 C (,M) '(*;') = p-tr\Xr l ) are the derivatives 

with respect to * of a(x) and (Q (rf -'>(jc) evaluated at x = , respectively. 

20 16. A method for correcting errors in data stored using Reed-Solomon encoding, the method 
comprising the steps of: 

a. using a dedicated circuit to compute syndrome values S,,S„. when tiu 

using 

forl<A<cz-l 

/•0 /«0 i-0 



25 »-| n-l n-l 

s.-z^-z^+z*** 



unless =0 for 

30 b "sing a dedicated circuit to determine a modified error locator polynomial 

p • o-(x) and a modified error evaluator polynomial fi ■ fi(x) from S k forlzkzd-l via an 

inversionless Berlekamp-Massey algorithm; 

c. using a dedicated circuit to compute the roots of /? • cr(x) using a Chien search, 
3 5 the roots of /? • cr(x) being inverse locations of the v errors; and 
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10 



15 



20 



25 



30 



35 



d. using a dedicated circuit to compute error values according to the equation: 
Y Q(*; ! ) D^ X \X] X ) • (50) 

17. A general purpose computer for correcting errors in data stored using Reed-Solomon 
encoding, the general purpose computer being configured to solve an error locator polynomial 
and an error locator polynomial by iteration without performing polynomial division and inverses 
during iteration. 

18. The general purpose computer as recited in Claim 17, wherein the general purpose 
computer is further configured to: 

a. define a modified error locator polynomial according to the equation: 

C ( '-'> (*) = p .C^Oc) = fi-a(x) = (fia Q9 fia l9 ...fia^) ; and (44) 

b. define a modified error evaluator polynomial according to the equation: 
D«-»{x) = p.D^\x) = f} •fiW = (/?.fi fl) ffi l ,..,/?.n,J. (45) 

19. The general purpose computer as recited in Claim 17, wherein the general purpose 
computer is configured to calculate the error locator polynomial and the error evaluator 
polynomial simultaneously without computing an inverse in a finite field. 

20. The general purpose computer as recited in Claim 17, wherein the general purpose 
computer is configured to iterate according to a modified Berlekamp-Massey algorithm. 

21. The general purpose computer as recited in Claim 17, wherein the general purpose 
computer is configured to: 

a. initially define: 

C i0) (x) = 1 >D {0) (x) = 1 ,Z'°> (x) = 1,5 (0) (jc) = 0, 7 (0) =0, 
=1 if Jt<0 ; 
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10 



15 



20 



25 



c. 



set k = k + 1 if k = d - 1 , otherwise, defining 



and let 

C<*)( X ) = r <*-» C<*"'>(x)- Z?*> Z<*-'>(x)x 
£ (i) (x) = r <*-»> B < *-»(x)-Zf*>ff< w »(x).xl4 



5 ( *>(x) = 
?<*>(x) = 



xy^x), if S« = 0 or if 2? ( *" I> > Jfc-I 

C ( *- ,} (x), if A (4) *0 and 2? ( *-'> s Ar-1 

^ ( *- l) (x), if A™ =0 or if 25 ( * l > > Jfe-1 

D<*' l \x) t if A ( *>*0 and <. k-l 



if A ( *> = 0 or if 2« ( *- |) > k-l 



? ( *" !) , 

.*-«<*-'>, .if A<*>*0 and 
<*- l > if A<*> = 0 or if 2? ( *- ,) >A:-l 



2? ( *-» s Jt-I 



(23) 



A w , if A ( *>*0 and 2{ ( *" ,) <. k-l 
; and 

d. repeat step (a). 
22. The general purpose computer as recited in Claim 17, wherein the general 
computer is further configured to: 

a. determine an error locator polynomial according to the equation: 
C<"-»(x) 



(24) 
(25) 

(26) 
(27) 

(28) 
(29) 

(30) 
(31) 

(32) 
(33) 



30 



<7(X)=- 



purpose 



(47) 



P 



;and 



determine an error evaluator polynomial according to the equation: 
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23. The general purpose computer as recited in Claim 17, wherein the general purpose 
computer is further configured to determine error magnitudes according to the equation: 

Y _ D id - l \XJ l ) (50) 

' no-*,*;') 

wherein ^ x - 1 ) = Xj Yl(l -X t XJ l ) and '(*;') = fia\X: x ) are the derivatives 
t*j 

with respect to * of a(x) and (Q (rf_,) (x) evaluated at * = , respectively. 

24. A general purpose computer for correcting errors in data stored using Reed-Solomon 
encoding, the general purpose computer being configured to: 

a. compute syndrome values 5", ,S 2 ^.S^ when teu using 

>«=0 /«=0 i>0 

unless 5 4 = 0 for l^*<rf-l; 

b. determine a modified error locator polynomial £ • a(x) and a modified error 
evaluator polynomial fid(x) from 5, forl£*;:</-l viaaninversionlessBerlekamp-Massey 
algorithm; 

c. compute the roots of fi ■ a(x) using a Chien search, the roots of p ■ a(x) being 
inverse locations of the v errors; and 



compute error values according to the equation: 
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Q(*;'> _ g(j£)_ d (</ - |) (^;') (50) 
n(l-Jf,*;') *;V(*;') *;'C<'-'>'(jr;') 



25. A device comprising a dedicated circuit configured to correct errors in data stored using 
Reed-Solomon encoding and further configured to solve an error locator polynomial and an error 
evaluator polynomial by iteration without performing polynomial division and inverses during 
iteration. 

26. The device as recited in Claim 25, further comprising a dedicated circuit configured to: 

a. define a modified error locator polynomial according to the equation: 

C <'-> (*) = 0 ■ C<"-'> ix)mft. cr{x) = (fia 0 , /?<r, ) ; and (44) 

b. define a modified error evaluator polynomial according to the equation: 
D iJ - l \x)^/}.D^\x) = p.n{x) = ^.Cl tt ,p.^..p.Cl J . l ). (45) 

27. The device as recited in Claim 25, further comprising a dedicated circuit configured to 
calculate the error locator polynomial and the error evaluator polynomial simultaneously without 
computing an inverse in a finite field. 

28. The device as recited in Claim 25, wherein the dedicated circuit thereof is configured to 
iterate according to a modified Berlekamp-Massey algorithm. 

29. The device as recited in Claim 25, further comprising a dedicated circuit configured to: 

a. initially define: 

C <°> (x) = 1 , 0 <°>0t) = I , A«\x) = 1 , F (0) (jc) = 0,7<°> = 0 , 
y (l) =1 if k<>0 ; 

b. set k = k + 1 if k = d - 1 , otherwise, defining 

a^zV-V ; (23) 



and let 
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30. 



C ,,, (x).y M C«-»(x)--8 l) A< k -»(x)-x 
> (x) = r (t - l) D' 1 -" (x) - A™I<*-*>(*) • x 14 



A«\x) = 
5<*>(x) = 

i <*>(*) = 



x-A ik -%), if S»=o or if 2i lk - l) >k-l 

C<*">(x), if A<*>*0 and 2!^ * k-i 

r5 (k - l \x), if A^O or if 2?**" 0 > k-\ 

D (k - l \x), if A ( *>*0 and 2« ( *' ,) s *-l 



? (*-t) 



if A w «0 or if 
if A (i) *0 and 



2? ( *" ,) s *-l 



(24) 
(25) 

(26) 
(27) 

(28) 
(29) 

(30) 
(31) 



Y<*> = 



Y ( *-'\ if 3»-0 or if 
A ( *>, if A w *0 and 



2t 



(*-i) 



> Ar-1 



(32) 
(33) 



2?<*- ,) s Jfc-1 
; and 

d. repeat step (a). 

The device as recited in Claim 25, further comprising a dedicated circuit configured to: 
a. determine an error locator polynomial according to the equation: 

C (-W) (x) (47) 



P 



; and 

b. determine an error evaluator polynomial according to the equation: 

/>(*) = n(x)-<r(x) = ^^-cr(x) * 

P 



(48) 



31. The device as recited in Claim 25, further comprising a dedicated circuit configured to 
determine error magnitudes according to the equation: 

J H<\-X t X?) X- l a\X- 1 ) X] l C«-»'(X- 1 ) 
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wherein &{X ?) = ^fjd " ) 3,141 C "->'(*;') = 0 ■ a\X^ ) are the derivatives 

with respect to * of a(x) and (Q^ l \x) evaluated at * = , respectively. 

32. A device for correcting errors in data stored using Reed-Solomon encoding, the device 
comprising a dedicated circuit configured to: 

a. compute syndrome values 5, , S 2 when tk u using 

/=0 J-0 /*>0 

unless =0 for I ^ Ar ^ ^ — 1 ; 

b. determine a modified error locator polynomial p a(x) and a modified error 
evaluator polynomial /} . Q(x) from 5, for 1 <: * <; rf - 1 via an inversionless Berlekamp-Massey 
algorithm; 

c. compute the roots of ft . a(x) using a Chien search, the roots of fi . <r(jc) being 
inverse locations of the v errors; and 

d. compute error values according to the equation: 

Y _ _C1(X: X ) Q(X?) D«~ l \x: 1 ) • (50) 



7 riO-Vj 1 ) X?<r\*l X ) X?C™'(X?) 



33. A device for correcting errors in data stored using Reed-Solomon encoding, the device 
comprising 

a. a first circuit performing iteration according to an inversionless Berlekamp- 
Massey algorithm; and 

b. a second circuit for defining a modified error locator polynomial according to the 
equation: 
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C<"-'>(x) = P C^\x) = fi-crix) = <Jkj„fkr x ,...p* d _ x y, and (44) 
for defining a modified error evaluator polynomial according to the equation: 

D lJ - l \x) = p.D^\x)^p.a(x) = {p.Cl,,p.n if ...p-Cl d . l ). (45) 

34. The device as recited in Claim 33, wherein the first and second circuits are defined by a 
very large scale integration (VLSI) chip. 
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